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\ ABSTRACT 

This thesis presents a method for estimating the distribution of 
stress components around the outer circumference of a cross section of 
a pipe, from the readings of strain gage elements arbitrarily positioned 
and oriented around this circumference. A least-squares procedure is 
used to obtain best estimates of the coefficients of Fourier expansions 
describing such distributions. A digital computer program was developed 
for applying the method. Data for testing the method and program were 
generated by a computer program Meine fhe best available theory of stress 
analysis in pipes. Methods for adding random errors to the data were 
adapted and used for closer simulation of actual situations. 

A second problem treated in this thesis is that of inferring the 
loading acting through the cross section of a straight pipe of concen- 
tric bore from known stresses at points on the external surface of the 
cross section which is presumed to be distant from stress concentra- 
tions. It is shown that this inference can be made from stresses at 


only three points of the cross section. A digital computer program was 


developed to do this. 
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Li. Introduction. 

In a thesis presented to the Naval Postgraduate School in May, 1966, 
LT Joseph W. Koch, Jr., USN, did pioneer work on the problem of obtain- 
ing optimal inferences of stresses in submarine main sea water piping, 
from strain gage readings taken during dockside hydrostatic tests and 
during submerged operations. This work provided a sequential procedure 
for determining such stresses, but it also indicated the need for many 
related investigations. 

Two of these are treated in the present thesis. The first is that 
of inferring, in an optimal way, the distributions of stresses around 
the outer circumference of a given cross section of pipe from observa- 
tions of strain obtained from an arbitrary distribution of strain gage 
elements around that circumference. There are several reasons for being 
interested in this problem. One reason is that even though the gage 
elements may be arranged in special patterns intended to simplify deter- 
mining stresses (for example three-gage rosette arrangements), failure 
of any of the elements might make it difficult to make effective use of 
valid information obtained from closely related gages. Another reason 
is the ability to infer stresses from arbitrarily arranged gages, which 
is a logical step in investigating the most effective and most economi- 
cal arrangements of gage elements. 

Briefly, this problem is solved herein by constructing Fourier 
series representing the desired stress components. Theoretical consider- 
ations are used to evaluate the coefficients of the series in an optimal 
way, minimizing the sum of the squares of the residuals. A digital com- 


puter program was developed to support the theory. 


The second problem treated here is that of inferring the force and 
moment components applied through a cross section of pipe and the inter- 
nal pressure acting inside the pipe. It is shown that these values may 
be inferred from a knowledge of the state of stress at only three abritrar- 
ily chosen points at the external circumference of the cross section, and 
a digital computer program for making the corresponding evaluations was 
developed. 

The solution of this problem assumes that the stress distribution is 
that given by what may be regarded as the best available current theory 
for the state of stress in a section of uniform pipe having concentric 
bore and remote from stress concentrations and other disturbances in 
the stress pattern. 

A digital computer program to generate numerical data from the best 
available theory was developed. This data was used to test and verify 
the theory and the computer programs. Suggestions of how to extend the 
present work are given in Section 4. 

The author would like to express his sincere appreciation to Dr. 
John E. Brock, Professor of Mechanical Engineering, Naval Postgraduate 
School, for his continued patience, efforts and most capable guidance 


while acting as faculty advisor in this work. 


Ze Least-Squares Estimation of Stress Distributions at the txternal 
Circumference of a Pipe by Means of Fourier Expansions, Using 
Strain Gage Data. 

The usual method of estimating the stress distribution at 
the surface of a body is by taking readings of strain gages grouped in 
rosettes, to determine local principal stresses. Unless the strain 
rosettes have a redundance of gage elements, when one or more of them 
fail, it is difficult to make effective use of the information provided 
by the elements which do not fail. 

The work that follows presents a method to estimate the stress 
distribution around an external circumference of a cross section of pipe, 
along a sequence of strain gage elements, which may or may not be group- 
ed in rosettes. If one or more gage elements fail, the information 
from the others is still useful. 


It is assumed that strain gage elements are arbitrarily placed 





(a) (b) 


Pig. 2.1. 


and arbitrarly oriented around the external circumference of a cross 


section of pipe so as to measure the strain produced by application of 


loading. Assume that there are N such elements as indicated in Fig. 
21a 
Define two orthogonal directions at the cross section, one normal 
to it, called longitudinal and represented by the letter L, and the 
other tangential to the external circumference, called tangential and 
represented by T; Fig. 2.1.b. Also define the angle 9 as the angle be- 
tween the vertical line that passes downward through the center of the 
cross section, and the radial line that passes through any specified 
point of the external circumference, measured as indicated in Fig. 2.1l.b. 
Measure the angular orientation of each gage element in the counter- 
clockwise direction from the tangential axis T, and call this angle @. 


A set of parameter-pairs can then be formed 


which describe the position and orientation of each element, Fig. 2.2. 





Bae. 2.2. 
At each one of the gage locations, there are tangential, longi- 
tudinal, and shear strains, ( € Ta? (<€ La? a i 
and (€ TL i respectively. The strains 
E . and € J along axis y and z, rotated y 
the angle @ from the tangential direction 


Fig. 2.3., are related to the local Pigs 2,00 


strains € 


aie L TL 
: € € les 
€. — cos @ sin @ 7 5; TL cos @ sin @ 
é S « e€ : 
pe ”" €. sin @ cos @ +5 TL 1} {| sin ) cos ¢ 


The corresponding relations for the strains ( € oa and ( € a 


along and transverse to the ith strain gage axis, are then 


Ck), = (<,),sin 6, + (€,), cos, f CO ean d,cos OF (2.1.1) 


2 2 
= © ; -(€ 
(é€ ) ( 4°08 a. + (€,),sin ¢g ( TL 18178 ,¢°5 0. @’, 1-2) 
At the external surface, strains are related to the corresponding 


stresses as follows 


¢. = (s,-v 8.) /E oil) 
€. =(s,-vs,) /E (2.2.2) 
€ 

T, = 8a iG (2.2.3) 


where s is stress 
v is Poisson's ratio 
E is Young's modulus of elasticity 
G is the shear modulus of elasticity 
Substituting these expressions for the strains in the relations 
(2.1) for the strains along and transverse to the ith gage element and 


simplifying 


Kew) = (3,,), [1 - (v + 1)sin’G, 1/8 + (s,), [1 - (v + L)cos“G, J/E 


+(s,,), sin 2 9. / (26) (273.1) 


(Gm) amanda - @v + L)cos 6, 1/E + (s,),[1 - Ww + 1)sin’9,1/E 
- (Spy), 8in2 8, /(26) (2.5352) 


A strain gage element is sensitive not only to the strain along 
its own axis but also, to a small extent, to the strain at right angles 
to this axis. The transverse sensitivity is usually hard to obtain. It 
is rarely available from the manufacturer and may have to be inferred from 
an experiment. Nevertheless, it may be of importance and it is included 
in the theoretical development that follows. 
For a particular gage element properly mounted on a given definite 
material, the relative change in electrical resistance, which is sensed 
by the strain indicator (or recorder), is given by a formula of the sort 
QOre=F € + F € (2.4) 
[ae yes, 
where R is the electrical resistance of the element, a and é. are 
the strains along and perpendicular to its axis, and Se and Fe are co- 
efficients which depend on the characteristics of the gage and the pro- 
perties of the material upon which it is mounted. 
The ratio 


k=F, / kes (2.5) 


is called the transverse sensitivity (factor) of the gage, and the 
quantity 


2 
F'=(1-k’)F 
( Dia 


is called the gage factor (which should be set on the sensor instrument). 
F' is the ratio of unit change in electric resistance to the quantity e, 


(indicated strain) sensed by the instrumentation. 
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Thus e= ( AR layer ' 


2 
er é, oie Ee /TQ - Esso) 
2 
= ( oy k€é, )/(1l - Kk’) (2.7) 


This is what the instrumentation should indicate if there were no 
errors. Clearly for a gage with no transverse sensitivity, k = 0, and 
e=€. 
y 
After substituting the expressions (2.3) for (és and ( € Oe in 


equation (2.7) and simplifying, the theoretical indicated strain e 


sensed along the ith strain gage element is given by 
e, = (s_), [1 - kv) - (L - (1 + v)sin 8, )/[EC1 - &7)) + 

+ (s,),[(1 - kv) - (1 - k)(1 + v)cosG, J/[E(1 - k°)2 + 

+ (8p,),(1 - K)sin2g,/[26(1 - k7)} (2.8) 


Let the coefficients of the stresses in the above expression be 


represented by 


AQ,) = A, = [(L - kv) - (1 - k)(1 4v)sin“6,1/[E - k*)] (2.9.1) 
B(O,) = B, = [(1 - kv) - (1 - k)(1 + v)cosG,1/lz (1 - ik“) 
(2.9.2) 
C(O.) = C, = (1 - k)sinag, /[2G(1 - k*)] (2.9.3) 
Then equation (2.8) can be put in the form 
e, = A, (85); + Bi (s,); + Ci (sp); (2.10) 


Any one of the stresses (Si) > (3,5 and (3.7) has continuous and 


finite values around the external circumference of the pipe cross 


ll 





stress 


Fig. 2.4. 


section, and can be represented as in Fig. 2.4. 


The functional relation between stress and the angular position 


@ can be expanded in an infinite Fourier 5 Oona 


(s,), = = ($4) 698 nd, + = (st) sin no, (2.11.1) 
- = =, 
(85); = (8,4) 698 no, + — (353), sin no, (2. Lies 
=, = 
(Sa) 5 = = (S74) 1°08 no, + = (8475) ,810 nO, (2.11.3) 
where 
(Sin » (8p), > (8p), > (Srp), >) GSouadn > 29d (Sopp), are 


as yet unknown numerical coefficients. 
Sufficiently accurate values can be obtained if the expansion is 
carried to a finite number of terms, say P terms, in which case the 


index n would reach a maximum value of P/2. Call this value M. Then 


” 4 
(3), =~ = (8.4) ,6°8 nd, + = (Srp) Sin n@, (2.12.1) 


The conventional expansion of a Fourier series is a constant term 
followed by a summation from 1 to o@ . The same is obtained by eliminat- 
ing the constant term and making summations from 0 to co . This latter 
form will be used in this work to facilitate mathematical manipulations 
and later programming for a digital computer. 


12 


Mz 


as 
=, (3,4) ,COs ne, OF 


(s1), ~ = o. ($73) ,8in 0, (2.12.2) 
s = 
(857) 5 ~~ = (Sara) nS nd. + oe (Sapp) Sin nd, (2.12.3) 


Note that this expansion contains 6(M + 1) Fourier coefficients. 
If these expressions for the stresses are substituted into equation 


(2.10) for the strain ess theoretically sensed by the ith strain gage 


element, 
= : 
—_—s 
ce A, | ba (S54) COS no, + — ree a 
m a | 
+B.{ Bo a ere) (2.13) 
1 W290 wre 
= 
+ c,[ = etd = (Sap) sin n@, 


Call ey the strain actually measured by the ith strain gage ele- 
ment, and E. the difference between that measured strain and eu» the 
strain that theoretically it should have measured, so that this error 


E. is expressed as 


E, =e, - e* (2.14) 


> 2 
_ oni (E,) 
a. Saenee2 
Ae fe, en) 
N 9 A ad 2 
— 9 @? - 2S ee ee cme) 255) 
421 i ae at esl 


To minimize R, partial derivatives respect to each of the 6(M + 1) 


Fourier coefficients are taken and set equal to zero: 


OR =0 
2D, 





where D, = 


c™ (ny > (pe > (rye > Cp), >» Spray > (Spy), in turn 


is 


Boe = Oh, Biss 4 iM 


Taking derivatives of R from expression (2.15) 

















v “»} e, N ne 
JR =0=2 BH e, = -2 SB e® ST 
dD, a - "O k ast rr) Dy 
Simplifying this 
N NV 
exe “Oo @. 
= = = ex : (2.17) 
4st ) k aS il 2 Ek 
Taking derivatives of e, from expression (2.13) 
goede = A, cos ke, fork =0, 1lje2,4.. ,M (2908e8) 
2 (Sp), 
de, 
=———~ = A,sin k@ (2.16.35 
2 (S58) i i z 
8 & 
ri = B cos ko. (2. 18.3) 
a %; 
> (Sip), ~ B,sin ke, (2.18245 
oe, 
= C.cos k@ (2.18.5) 
Ce) i . 
TLA’ k 
8 & 
> (Spr5) = C,sin ke, (2.18.6) 


Note that there are 6(M + 1) derivatives of this form. Substitut- 


ing in equation (2.17) the expression for e, from equation (2.13) 


Zf, p= (s,, LB) net no, +..f+..+¢,[ 


oe; 2 &% 
Ace = (Srp) Sin lies = = ex 2D, 
(2.19) 
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which is the same as 


M ¥ ae w 1. oe Bee 
Bou, Z Aone 28] + BGu), Baar 3S 1 


| 4 ee | 7 te oe. 
i ae I Cc" le 
+- =e o +> = (Srza) 2 0 AHA n&. oy “ce < > D (2.20) 


pee ke=s0H lyeZ2, «=. , M 
When the index n in the Fourier expansions is equal to zero, the 


terms containing sin(n6, ) are zero. This makes the coefficients (S55)? 


(s13) and (8.73) also zero, reducing the Fourier i to (6M + 
3) and the same number of derivatives of the form os - There will 
k 


be (6M + 3) equations in-seach of which appear all the Fourier coeffici- 
ents. This makes a set of (6M + 3) equations with the same number of 


unknowns, (the Fourier coefficients), which can be solved. 


Once the Fourier coefficients are solved for, they can be substitut- 
ed in expressions (2.7) for the stresses. The stress distributions Sips 
sy and Say, are defined then as functions of the angle 9 over the com- 
plete external circumference of the cross section. 

A digital computer program was developed to apply the theory of the 
preceding text. It was tested using a mathematical model of a canti- 
levered pipe. The program is called TERESITA and is presented in Ap- 
pendix B of this thesis. 

It was assumed that strain gages having a zero transverse sensiti- 


vity factor were placed at the external circumference of cross section 


No. 1, see Fig. A.1. in Appendix A. Four arrangements of gage elements 


YD 


and several numbers of Fourier coefficients were used. 

The first arrangement consisted of twelve strain gage rosettes 
equally spaced 30°from each other around the cross section. Each rosette 
has three gage elements oriented at 0$ 45°and 90°from the positive tangen- 
tial axis. Runs were made with 11, 9, 7, 5 and 3 Fourier coefficients 
for each stress component. 

The second arrangement consisted of twelve rosettes spaced as above 
but with elements oriented at 0°, 60°, and 120°. Runs were made for the 
same number of Fourier coefficients. 

The third arrangement consisted of eight strain rosettes equally 
spaced 45° from each other. The rosettes had four gage elements orient- 
ed at 0°, 45°, 90° and 135° from the positive tangential axis. Runs were 
made for the same number of Fourier coefficients. 

As input data for testing the program, stress distributions pro- 
duced by the loading were evaluated at the external circumference of 
cross section No. 1. From the stresses, the theoretical strain that the 
gage element would sense, was computed. The stress distributions re- 
covered by the program agreed in all cases with the input information. 
Results using eleven Fourier coefficients (for each stress component Bs 


s. and 81) are almost identical with those using only three coefficients. 


L 
Additional runs were made using the second arrangement of gages, 
but eliminating arbitrarily the readings of some of them. This is what 
would have to be done in an actual case if some gage elements were known 
to have failed. Results are satisfactorily close to the preceding ones. 
Other runs were made using a strain gage transverse sensitivity factor 


of 1%. The program recovered likewise the input information. 


The fourth gage arrangement had one gage element every 10° around 
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the circumference with different orientation for each gage. Kung were 
made using the same number of coefficients. 

All the runs described above were repeated introducing in the gage 
element readings small random errors so to simulate actual situations. 
Calculated loadings differed from the input loadings to a degree which 
depended in each case on the amount of error introduced. 

The program and some sample numerical results of the runs made are 
presented in Appendix B of this thesis. 

A more detailed discussion of the effects of changing the number 
of Fourier coefficients and gage element arrangement is presented in 


Appendix E. 
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Di Inference of the Loading Acting at a Cross Section of a Pipe of 
Concentric Bore from Stresses at the External Circumference of 
that Cross Section. 

If stresses are known at a sufficient number of points around 
the external circumference of the cross section of a pipe of concentric 
bore at a distance from stress risers and concentrated loads, the load- 
ing acting at the cross section can be determined from theory. 

It is intended to analyse the effects that such loads produce 
at a general point P under the conditions specified above. From this 
analysis, general conclusions will be reached which lead to a mathe- 
matical procedure to determine the loads. 

Consider a cross section of pipe as specified above, Fig. 3.1. 


Choose reference axes x, y and z with the center point 0 as origin. Let 










ee 


eat | 






(a) (b) 


Pape oe 


i, J and k be unit vectors along the positive x, y and z axis so to form 


a right handed triad. Let eR be a unit vector along the outwards direc- 


tion of the radius that passes through the point P, and en another unit 


vector so that & and k form a right handed triad. The angle @ is 


R’? “T 


measured from i counterclockwise to ioe Consider the differential area 


surrounding the point P. Assume that tensile stresses Sp and Sy (along 


the directions of 6, and k respectively) and shear stress s 


duced by internal pressure p and by a general force F and a general mom- 


are pro- 
tH i 
ent M acting on the cross section at the center point 0. 


Let the components of F and M in the @., @. and k triad be (FR> F 


R’ T T 
F) and (M, » Mp M,) respectively, and in the i, j, k triad, (FL co a 
and (Ms MoE M)- 


From strength of materials considerations reviewed in reference 5, 


the stresses produced by p, F and M at the external point P are: 


Te 2p ~ pL Se - a) (Seb) 

S., = -M,, b/It FS JA (322) 

Sart, = M, b / (2 I) +D FS (3.3) 
where a = internal radius of the cross section 


b = external radius of the cross section 
A = cross sectional area 


I = moment of inertia of the cross sectional 
area about a diameter 


Tey RG 2) on ee 


v = Poisson's ratio 


n =v + 3/2 


Making use of the two dimensional coordinate transformation tensor 


M, cos 0 sin 0 M 
x 

M -sin 9 cos Q M 
T y 


M. can be expressed as 


Mh = - M sin 9 sary cos 9 (3.4) 
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Similarly, 


Fe = aha sin @ + Fy cos Q 325) 


The expression (3.2) for 8, can be transformed to 
he (1/A)F | + (b/I)sin 0 M - (b/I)cos 0 M (3.6) 
Multiplying through by (I/b) 
(I/b)s, = (I/(Ab))F, + (sin Q)M, - (cos QM, (3.7) 


Letting (I/b) = g, and (I/(Ab)) = h, both of which are cross 


sectional constants, substituting above and changing sides 
+ - - 
h +. sin @ M cos Q M.. gs, (3.8) 


The preceding theory is applicable to a point on the exterior 
surface of a piece of straight, uniform, concentric pipe at a suf- 
ficiently great distance from its ends or other connections which might 
introduce local disturbances into the stress distribution. Presuming 
that these conditions are satisfied, we next turn to the matter of try- 


ing to infer the loading conditions at the cross section, i.e., Pye Ea 


Foe Mi» Ms M, and p, from a knowledge of the stresses at several points 


on the external surface of the cross section. 


For a general point P, s, and 90 are known and the only unknowns in 


L 
equation (3.8) are Fo M and ay! Taking three such points at the cross 


section, the three unknowns can be determined proceeding as follows. 


1’ 8. and 
8. and the axial stresses are 81> 8, and 84 respectively. A set of 


equations can then be formed 


Assume that the angular positions for the points are 9 


h F. + s® M - c8, ~ = gs (3.9.1) 


1 1 


h F + sQ M. - cQ 7 = £8, (3.9.2) 


2 2 
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h Be + sQ M, - c@., i = 9s (3.9.3) 


3 3 


where the letters s and c with no subscripts are abreviations for sin 
and cos respectively. These abbreviations will be used throughout this 
section. 


The set of equations (3.9) is solved applying Cramer's rule 


h $8, -cO, 
Denominator = h 88. -c@, 
h 38., -c@., 


Bs) sO -cQ 


1 1 
Numerator of F, * 8s. SQ, -c8, 
gs. s®., ~c®., 


= -2( $,5(8, - 8.) + s,5(0, > 81) a s,5(8, rm 8.)) 


(3.11) 
h gs) ~c8, 
Numerator of M #= h gs -cQ 
x 2 2 
| h gs. -c®, 


= -hg( 0, (8, - 8.) + c8,(8, E s.) + c8,.(s, - 8,)) 


(3.12) 
h s8, gs, 
Numerator of My, = h 58, gs. 
h s0, gs, 


= -hg (80, (8, = $5) + s0,(8, 2 35) + s0,(s. - s,)) 


(3.13) 
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The values of F M, and bd are found by dividing the respective 
numerator by the denominator. 

The remaining components of the vectors F and M appear in the 
expression (3.3) for Snr? which by using equation (3.5) can be trans- 
formed to 


s = (b/2I)M_ - Dsin © F, + Deos 0 K (3.14) 


TL 
Multiplying through by (1/D), letting (-1/D) = g*, (-b/2ID) = h*, 
both cross sectional constants, substituting in (3.14) and changing 
sides 


h* M + sin®F - cos @F_= = gxg (3.15) 
Z x y 


TL 
This expression is of the same form as (3.8) which led to the 
evaluation of Fo M,. and Be by using information from three points. In 

this case the unknowns are M,> F and ey which can be evaluated by 
using the same three points and applying the same method. In this form 
the loads F and M are completely determined. 

Since 8. is produced by the internal pressure p alone, it can be 


solved for from equation (3.1) with the information from any one of 


the three points used to evaluate F and M 


p = ,(b~ - a°)/(2a°) (3.16) 

From the above development, it can be concluded that theoretically 
the loads F, M and p acting at the cross section can be determined from 
known stresses at three points of the external surface of the cross 
section. 

A digital computer program called TIBISAY was developed to apply 
the preceding treatment. The program was tested using a mathematical 
model of a cantilevered pipe as in Section 2 of this thesis. Stresses 


generated at cross section No. 1, at points 0°, 120° and 240° from the 
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vertical downwards directicn were used. The recovered loads were almost 
identical with the input ones producing the stresses. Deviations are 
due to small round-off errors introduced by the computer. The program 


and results of the test are presented in Appendix C of this thesis. 
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4. Conclusions. 

The results obtained in the tests for applicability of the theories 
developed in this thesis indicate that the proposed methods can be used 
for analysing stresses in pipes. 

There are, however, many factors influencing the validity of the 
results. Some of them are related to the amount of error introduced 
when determining geometries and material properties. Others are due to 
proximity to concentrated loads and stress concentrations. Other fac- 
tors are directly related to the degree of redundance of strain gage 
data, the number of gages used and their arrangement around the cross 
section, and the number of Fourier coefficients to be used. Also, it may 
be that the least-squares method used does not produce the best results; 
perhaps, by minimizing some other function of error, better results could 
be obtained. 

First steps were taken to investigate further some of these factors. 
Appendix E of this thesis presents results of a survey made concerning 
the influence of the number of Fourier coefficients used in the accuracy 
of results. There are also preliminary results of the influence of the 
number of gage elements and their arrangement around the cross section. 
Due to time limitations, further investigations were not conducted. 

The digital computer programs included in this thesis can be used 
to investigate deeper into the field. Some refinements of the methods 
can be added. It is believed that the following points should not be 
neglected. : 

a. Other criteria for minimizing errors can be investigated. 
A discussion of such criteria is given in reference 7 of this thesis. 
It is likely that some of the methods described there may have applica- 


tion in the pipe stress problem considered herein. 
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b. One may investigate the question of the minimum number of 
strain gage elements necessary to obtain acceptable stress distributions, 
and their orientation around the cross section. Similarly, the optimum 
number of elements and preferred orientation could be investigated. 

c. While observing the results for a determined gage arrange- 
ment, it was noticed that good results were obtained using from 3 to 7 
Fourier coefficients, but more than that produced bad results; using 
again 5 coefficients and changing the gage arrangement, good results were 
obtained for some arrangements and bad results for others. Some of these 
results are presented in Appendix G of this thesis. Further study could 
be made of the optimal number of Fourier coefficients to be employed. 

d. A method for data rejection should be incorporated in the 
theory and computer programs. There are instances when it is obvious 
that a strain gage element has failed and the reading of that gage should 
be neglected. Usually bad gages are not easily identified and bad data 
is processed along with good data. There is the possibility of incorporat- 
ing a data rejection procedure in the method developed in Section 2 of 
this thesis. One suggestion for such a procedure is given in the follow- 
ing. 

The estimated Fourier coefficients of the stress distributions are 
evaluated from a minimization of the sum of the squares of the errors, 
i.e., the differences between the gage element readings and what they 
should have read theoretically. We can assume that these errors are 
random having a statistical normal distribution with mean and standard 
deviation which could be computed. Establishing a level of significance 
for a deviation to be sufficiently apart from the mean, one can reject 
the data having greater deviations and recompute the stress distribu- 


tions with the remaining data. Similar screenings can be done with the 
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newly computed stress distributions until the remaining data have de- 
viations under the level of rejection. This theory could be readily 
used if the redundance of data were large. Usually this is not the case 
and one must avoid rejecting data which are not perfect but which are 
useful and necessary. If the data is known to have a normal curve with 
a very sharp peak, bad data is easily determined and low significance 
levels can be used. If the bell is fairly flat, the majority of the 
deviations have pronounced deviations from the mean and the level of 
rejection becomes difficult to establish. The writer is inclined to 
think that there are criteria that may satisfy all conditions. Once 
they are found, the theory could be applied to this problem and to 


similar ones. 
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APPENDIX A 
DATA GENERATION THEORY AND RELATED DIGITAL COMPUTER PROGRAMS 
The data used to test the applicability of the theories developed 
in this thesis was generated for a mathematical model of a cantilevered 
pipe of concentric bore subjected to internal pressure p and general 
force F and general moment M applied to the tip. Details of this pipe 
are shown in Fig. A.1. The following material properties, pipe dimensions, 


and loads were specified: 


i 





Fig. A.l. 
E = 30,000 ksi 
v = 0.3 
Specific weight = 0.283 iho see 


External diameter = 12.0 in 


Internal diameter 10.0 in 


F =F =F _ = 100.0 1b 
Xx y Z 
M =M =M _ #= 100.0 1b-in 
X y Zz 
p = 100.0 psi 
Subroutine CANTIL presented herein in this Appendix, was used 


to evaluate the tangential, longitudinal and shear stress distributions 
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at the external circumference of cross sections 1 and 2 of this pipe. 
The stress distributions were used in turn to supply the data required 
for the tests specified in other parts of this thesis. 

This subroutine is used by either of the (main) programs TERESITA 
or TIBISAY described in Appendices B and C. The arrangement of the card 
input is indicated in Fig. B.1. of this thesis. Description and specifi- 


cations for CANTIL follow this page. 
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DESCRIPTION OF SUBROUTINE CANTIL 


Subroutine CANTIL evaluates the state of stress at any point of 


a cantilevered straight pipe of concentric bore, with internal pressure 


and loaded at the tip by general forces and general moments. The program 


is capable of analysing from 1 to 360 points in each of 1 to 5 cross 


sections of the pipe. The points have to be at the same radial distance 


from the center of the pipe, and separated from one another by.angular 


intervals which are submultiples of 360°, from 1° to 360°. 


The data deck,see Fig. B.1,includes the following cards: 


Card No. 


a. 


Card No. 


1. Format 3F20.8 

External diameter of the pipe, inches. 

Internal diameter of the pipe, inches. 

Pipe length, inches. 

2. Format 3E20.8, 2F10.8 

Young's modulus of elasticity, psi. 

Shear modulus of elasticity, psi. This can be 
omitted in which case it will be computed by the 
program from E/2( 1+ v). 

Coefficient of thermal expansion. This can be omitted. 

It is not actually used in any present application. 

Poisson's ratio. 

Specific weight of the material, pounds per cubic 

inch. 

3. Format 215, 5F10.5 

Number of points to be analysed in each cross section. 

Number of cross sections to be analysed in the pipe. 

The distances from the tip of the pipe to the cross 


sections to be analysed, inches. 
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Card No. 


a. 


Card No. 


a. 


Card No. 


a. 


Card No. 


a. 


Card No. 


a. 


4. Format F80.2 
Radial distance of the points to be analysed, inches. 
NOTE: for applications to other parts of this thesis, 
evaluations are made only at the outside sur- 
face, but the program is capable of making 
evaluations at interior points. 
5. Format F80.2 
Internal pressure, psi. 
6. Format 3F20.8 
Components of the force applied at the tip of the 
pipe, vertical, sidewards and outwards from the 
tip, so as to form a right handed set of orthogonal 
axes, lbs. 
7. Format 3F20.8 
Components of the moment applied at the tip of the pipe, 
in the same arrangement as for the forces, lb-inches. 
8. Format F80.2 
This is a sentinel card with one number onit. A 
negative number indicates that the job is finished; 
a zero or any positive number indicates that the 
program will be rerun for another set of radial 
distance and loads, in which case, cards similar to 


No. 4 through 8 must follow this one. 


OUTPUTS (all arguments are outputs) 


RO. External radius of the pipe. 


RI. Internal radius of the pipe. 


XE. Young's modulus of elasticity. 
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XL1. 


NPOINT. 


AREA. 


ALL1. 


Shear modulus of elasticity. 


Poisson's ratio. 


Specific weight of the material. 


A constant evaluated when analysing shear stresses pro- 


duced by shearing forces. 


A one dimensional array containing the distances from 


the tip to each cross section considered. 


Number of points analysed in each cross section. 


Cross sectional area of the pipe. 


A three dimensional array intended to store the stresses 


at any of the points, in case the points are at the 


external surface of the pipe; otherwise, it has no mean- 


ing. 


Other subroutines which are needed and listings for which are given in 


this Appendix are: 


TORS : 


BEND: 


AXTAL: 


SHEAR: 


CROSS: 


DOT: 


TIMES ;: 


PLUS: 


MINUS: 


TRANSF: 


used by CANTIL to compute shear stresses produced by 


torsional moments. 


used by CANTIL to compute bending stresses. 


used by CANTIL to compute stresses due to axial forces. 


used by CANTIL to compute shear stresses produced by 


shear forces. 


computes 
computes 
computes 
computes 
computes 


computes 


the 


the 


the 


the 


the 


the 


cross product of two vectors. 

dot product of two vectors. 
product of a scalar and a vector. 
sum of two vectors. 

difference between two vectors. 


components of a vector in another angularly 


displaced set of axes. 
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METHOD. 

The vectorial approach discussed in reference 5 of this thesis is 
used to compute stresses produced at points of the cross section by 
internal pressure and forces and moments acting at the centroid of that 


cross section. 


M 





Fig. A-2. 
Those forces and moments, in turn, are computed by the classical vectori- 
al analysis of the statics of a beam. If the force F and the moment M 
are applied at one end of a beam, Fig. A.2.a., the force ie and moment 
M, acting at the centroid of the cross section nn, Fig. A.2.c., are 
found by considering the equilibrium of the segment of beam from the 
tip to the cross section. Summing forces and taking moments about the 


centroid of the cross section 
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Bacar = 0 (A.1) 


M + M + LxF + 43LxW = 0 (A.2) 


where L is the vector distance from the centroid of the cross section 


to the point of app 


lication of the loads at the tip, and W is the force 


due to the weight of the segment of beam. 


Solving for Be 


and M 

n 
oe = - (F + W) (A.3) 
Me = - (M + LxF + 4LxW) (A.4) 


The forces and moments acting at the cross section under considera- 


tion on the remaini 


ng portion of the beam, are the negatives of Ea and 


M, found above, Fig. A.2.b. 


PRINT -OUTS 


The print-outs are several pages as follows 


lst page. 
a. 


b. 


2nd page. 


Pipe dimensions 

Computed cross sectional area, moment of inertia 

and polar moment of inertia of the cross section. 
Material properties, 

Number of points analysed in each cross section, 
number of cross sections analysed and their distances 


from the tip. 


Radial distance of the points analysed . 
Loads applied to the tip. 
Distance from the tip to the first cross section 


analysed. 
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d. Computed loads acting at that cross section 

Following pages. For each analysed point the print-out contains: 

a. Angular position of the point in the cross section 
(degrees) 
b. Forces and moments acting at that cross section in 
radial, tangential and axial components. 
ec. Stresses produced by individual loads. 
SHTZT is the shear stress in the ZT direction 
produced by the axial component of the moment 
SBZZ axial (ZZ direction) stress produced by 
the tangential component of the moment. 
SAZZ axial stress produced by the axial compon- 
ent of the force. 
SHFZR and SHFZT are stresses produced by the 
radial and tangential components of the force. 
SPRR and SPTT are radial and tangential stresses 
produced by the internal pressure. 
ds If the analysed points are in the external surface 
of the pipe, there will be additional print-outs of 
the total axial, tangential and shear stresses at 
the point. 

Analysis of other cross sections for the initial and/or other radial 
distance and tip loads and pressure produce similar outputs as described 
from page 2 on. The complete program and sample outputs of the analysis 
made for points at the external surface and 4% in. under that surface of 


cross sections 1 and 2 of the pipe shown in Fig. A.1., follow this page. 
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APPENDIX B 

DIGITAL COMPUTER PROGRAM FOR SECTION 2 AND NUMERICAL RESULTS 

The programs presented herein provide a method for applying the 
theory developed in Section 2 of this thesis. Actual situations can be 
handled by subroutines ZULIA, GUAYRA, CUMANA and DSIMQ (DSIMQ is present- 
ed in Appendix D). Simulated situations require in addition subroutines 
CANTIL presented in Appendix A, and TIUNA and BLANCA presented in this 
Appendix. The last of the subroutines mentioned is an adaptation of a 
subroutine taken from reference 7 of this thesis. 

ZULIA is used to specify the required information, such as material 
properties, number of Fourier coefficients to be computed and descrip- 
tion of the arrangement of strain gage elements and their transverse 
sensitivity factor. Enough self explanatory comment statements have been 
included which specify the introduction of the information. Print-out 
statements present the information received for processing. Once the in- 
formation is properly arranged, subroutine GUAYRA is used to make the re- 
quired computations. 

The sequence in GUAYRA follows closely the theoretical development 
of the method. It has sufficient comment statements separating blocks 
of computer processing instructions as to follow the theoretical pro- 
cedure. Print-out statements present the computed Fourier coefficients. 

CUMANA is used by GUAYRA for auxiliary computations. 

Use of the subroutines is made by a very brief main program TERESITA. 
This program merely dimensions the various subroutines and turns control 
over to subroutine ZULIA. The card arrangement of the source deck is as 
shown in Fig. B.1l. 

Some sample results of the application of the program as specified 


at the end of Section 2, follow the card listings of the programs. 
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APPENDIX C 

DIGITAL COMPUTER PROGRAM FOR SECTION 3 AND NUMERICAL RESULTS 

Subroutine MACUTO presented in this Appendix provides a method for 
applying the theory developed in Section 3. Data for simulated situa- 
tions can be obtained from subroutine CANTIL presented in Appendix A. 

MACUTO has been programmed following closely the development of the 
method. Self explanatory comment statements are included, which detail 
the type of data to be used and the procedural sequence of the method. 
Print-outs present the data processed and the computed loading. 

Use of the subroutines is made by a brief main program called 
TIBISAY. This program merely dimensions the subroutines, accepts and 
prepares the input data, and turns control over to subroutine MACUTO. 
The card arrangement of the source deck is similar to that shown in 
Fige B.1. 

Some sample results of the application made of the programs as 
specified in Section 3 are presented following the card listings of the 


programs. 
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APPENDIX D 


Subroutine DSIIQ of the Comnuter Facility Center of the Naval Fostgrraduate 
-chool, to solve linear-simultaneous equations. 
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APPENDIX E 

Effects of Changing the Number of Fourier Coefficients and the 
Amount and Quality of the Data in the Application of the Theory Developed 
in Section 2. 

The theory developed in Section 2 of this thesis estimates Fourier 
coefficients of the distributions of axial, tangential, and shear stress- 
es assuming that enough data has been obtained by placing a sufficient 
number of strain gage elements around the external surface of a cross 
section of pipe. The correctness of the estimated coefficients depends 
on their number and on the amount and quality of the data. Depending on 
the relative orientation of the gage elements, they can pick information 
useful for determining the coefficients of all the distributions, or, if 
they are arranged in certain disadvantageous ways, they may be able to 
provide information for none, or oniy some of the coefficients. 

The requirement of inferring a given number N of Fourier coeffi- 
cients, demands having, at a minimum, a corresponding number M = M(N) 
of valid data items. If more data are available, there is a redundance 
of data which permits more accurate determination of the N coefficients. 

Specific rules for the necessary number and disposition of gage ele- 
ments to obtain a specified number of correct coefficients are not avail- 
able, but they probably could be worked out. Some idea of the gage ele- 
ments requirements can be obtained by testing actual cases and comparing 
results with correct ones obtained from other methods. The digital com- 
puter programs included in this thesis can be used as they are to perform 
such tests. The procedure tc follow .is to simply specify in subroutine 
ZULIA the number of gage elements, their angular positions and angular 
Orientations around the cross section, and the number of Fourier co- 


efficients to be computed. The print-out of this program presents the 
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distributions obtained in this way, shown adjacent to the theoretically 


correct ones, so that comparisons can be readily done. 


A preliminary survey of gage elements arrangements was made using 


that procedure, applied to data with no added random errors. Results 


obtained are presented following these pages. Many of these results 


might have been anticipated from theoretical considerations. The follow- 


ing nomenclature is used when applicable: 


(a) 


(b) 


No. of coefficients 


3 to ll: the results presented were obtained using from 
three to eleven Fourier coefficients per series. These 
results apply to that number of coefficients. 

73; results apply only to seven coefficients. 

Results 

good: indicates that the coefficients obtained describe 
correctly the tangential, axial and shear stress distri- 
butions. 

bad: The computed coefficients do not describe either 
of the distributions. 

only shear good: good coefficients were obtained only for 


the shear stress distribution. 


Inasmuch as the data were without added errors, the recovery should be 


and actually was excellent or no good at all; accordingly, there is 


nothing subjective in the judgment "good" or "bad". 
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Gage elements located at 12 points 
equally spaced 30° around the external 


circumference of the cross section. 


gage elements 


arrangement No. of Fourier 
at each point coefficients per series results 
45° 
Ma ee 3 to ll good 
T axis 
90° 
oe: 3 to li good 


N59 3 to Ll only shear good 
Ais 3 to Li bad 


ay £0211 bad 


3 to ll good 


60° 
\ fo 3: £Oui: only shear good 
foo" 3 to ll bad 


re 


II) 
Gage elements located at 8 points 
equally spaced 45° around the external 


circumference of the cross section. 


gage elements 


arrangement No. of Fourier 
at_each point coefficients per series results 
45° 45° 
\ bis" SCOny good 
9 and 11 bad 
same 
90° 


Yo: 3 tom? good 


Samia 9 and lil bad 


90 
Yo: 3 46Gr, only shear good 


same 9 and ll bad 
f 45° 3 to tl bad 

90° 
3 tock bad 
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IIT) 


IV) 


One gage element every 10° around the external 
circumference of the cross section, totaling 
36 elements. The first element at 0° from the 
tangential direction, and every consecutive 


one, incrementing this angle by 10°. 


from 3 to 11 coefficients: bad results 


One gage element every 10° around the external 
circumference of the cross section, totaling 
36 elements. The first element at 0° from the 
tangential direction, and every consecutive 


one, incrementing this angle by 50°. 


from 3 to 9 coefficients: good results 


11 coefficients: bad results 
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